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In the context of f(R) theories of gravity, we study the evolution of scalar cosmological perturba- 
tions in the metric formalism. Using a completely general procedure, we find the exact fourth-order 
differential equation for the matter density perturbations in the longitudinal gauge. In the case of 
sub-Hubble modes, the expression reduces to a second-order equation which is compared with the 
standard (quasi-static) equation used in the literature. We show that for general f(R) functions 
the quasi-static approximation is not justified. However, for those functions adequately describing 
the present phase of accelerated expansion and satisfying local gravity tests, it provides a correct 
description for the evolution of perturbations. 
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INTRODUCTION 

The present phase of accelerated expansion of the uni- 
verse [l| poses one of the most important problems of 
modern cosmology. It is well known that ordinary Ein- 
stein's equations in either a matter or radiation domi- 
nated universe give rise to decelerated periods of expan- 
sion. In order to have acceleration, the total energy- 
momentum tensor appearing on the right hand side of 
the equations should be dominated at late times by a 
hypothetical negative pressure fluid usually called dark 
energy (see Q and references therein). 

However, there are other possibilities to generate a pe- 
riod of acceleration in which no new sources are included 
on the r.h.s. of the equations, but instead Einstein's 
gravity itself is modified Q. In one of such possibili- 
ties, new functions of the curvature scalar (f{R) terms) 
are included in the gravitational action, which amounts 
to modifiying the l.h.s of the equations of motion. Al- 
though such theories are able to describe the accelerated 
expansion on cosmological scales correctly, they typically 
give rise to strong effects on smaller scales. In any case 
viable models can be constructed to be compatible with 
local gravity tests and other cosmological constraints 0] ■ 

The important question that arises is therefore how to 
discriminate dark energy models from modified gravities 
using present or future observations. It is known that by 
choosing particular f(R) functions, one can mimic any 
background evolution (expansion history), and in par- 
ticular that of ACDM. Accordingly, the exclusive use 
of observations such as high-redshift Hubble diagrams 
from SNIa l|, baryon acoustic oscillations [B| or CMB 
shift factor [6| , based on different distance measurements 
which are sensitive only to the expansion history,cannot 
settle the question of the nature of dark energy [7J] . 



However, there exists a different type of observations 
which are sensitive, not only to the expansion history, 
but also to the evolution of matter density perturbations. 
The fact that the evolution of perturbations depends on 
the specific gravity model, i.e. it differs in general from 
that of Einstein's gravity even though the background 
evolution is the same, means that this kind of observa- 
tions will help distinguishing between different models 
for acceleration. 

In this work we study the problem of determining 
the exact equation for the evolution of matter density 
perturbations for arbitrary f(R) theories. Such prob- 
lem had been previously considered in the literature 
(EI EM , EEl) an d approximated equations have 
been widely used. They are typically based on the so 
called quasi-static approximation in which all the time 
derivative terms for the gravitational potentials are dis- 
carded, and only those including density perturbations 
are kept From our exact result, we will be able to 

determine under which conditions such an approximation 
can be justified. 

The paper is organized as follows: in Section 2, we 
briefly review the perturbations equations for the stan- 
dard ACDM model. In Section 3 we obtain the perturbed 
equations for general f(R) theories. In Section 4 we de- 
scribe the procedure to obtain the general equation for 
the density perturbation. In Section 5 we summarize the 
main viability condition for f(R) theories. Section 6 is 
devoted to the study of the validity of the quasi-static ap- 
proximation. In Section 7 we apply our results to some 
particular models and finally in Section 8 we include the 
main conclusions. In Appendices I and II we have also in- 
cluded complete expressions for the relevant coefficients 
of the perturbation equation. 



DENSITY PERTURBATIONS IN ACDM 
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Let us start by considering the simplest model for dark 
energy described by a cosmological constant A. The cor- 
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responding Einstein's equations read: 



PERTURBATIONS IN f(R) THEORIES 



(1) 



where G^ v is the Einstein's tensor and T^ u is the energy- 
momentum tensor for matter. 

In the metric formalism for the ACDM model it is pos- 
sible to obtain a second order differential equation for the 
growth of matter density perturbation 5 = Sp/po- Let 
us consider the scalar perturbations of a flat FRW metric 
in the longitudinal gauge and in conformal time: 

ds 2 = a 2 {r])[{l + 2&)drf - (1 - 2^)(dr 2 + r 2 dil\)] (2) 

where $ = $(r?,"a?) and "J" = ^(77, x ) are the scalar 
perturbations. From this metric, we obtain the first-order 
perturbed Einstein's equation: 



SG"„ = -87rGtfT"„ 



(3) 



where the perturbed energy-momentum tensor reads: 

<5T° = S P = Po S , ST^ = -5PS l 3 = -c 2 s 5 l ]Po 5 

ST , = -6T\ =-(l + c 2 s )p d. l v (4) 

with po the unperturbed energy density and v the po- 
tential for velocity perturbations. We assume that the 
perturbed and unperturbed matter have the same equa- 
tion of state, i.e. SP/6p = c 2 = Po/po, where c s = for 
matter perturbations. The resulting differential equation 
for 5 in Fourier space is written as: 



6" +H 



k 4 - 6pk 2 - 18p 2 



6' 



k 4 - p(3k 2 + 9H 2 ) 
„k 4 + 9p(2p - m 2 ) - k 2 (9p - 3H 2 ) 



-P- 



m 2 



6 = (5) 



k 4 - p(ik 2 

where p = AnGpoa 2 = —H' + 7i 2 and H = a' /a with 
prime denoting derivative with respect to time rj. We 
point out that it is not necessary to explicitly calculate 
potentials $ and \& to obtain equation but algebraic 
manipulations in the field equations are enough to get 
this result. In the extreme sub-Hubble limit, i.e. kr\ 1 
or equivalently k ^> H, ©is reduced to the well-known 
expression: 



5" + HS' - AnGp a 2 6 = 



(6) 



In this regime and at early times, the matter energy 
density dominates over the cosmological constant and it 
is easy to show that 5 solutions for © grow as a(rf). At 
late times (near today) the cosmological constant contri- 
bution is not negligible and power-law solutions for ([6]) 
no longer exist. It is necessary in this case to assume 
an ansatz for S. One which works very well is the one 
proposed in and [HI 



5(a) _ J2 [n m (a)^-l]dlna 

a 



(7) 



This expression fits with high precission the numerical 
solution for S with a constant parameter 7 = 6/11. 



Let us consider the modified gravitational action: 
1 



S 



16?rG 



d 4 Xy/=g~(R + f(R)) 



(8) 



where R is the scalar curvature 1 . The corresponding 
equations of motion read: 

Gfj,u - ^g^f(R) + R^f R (R)~g^Df R (R) + f R (R), tw 



(9) 



where f R (R) = df(R)/dR. For the background flat 
Robertson- Walker metric they read: 



and 

\(H' + 2H 2 )(l + / fl )-i(i? + /o) 
a z 2 



-87rGpo (10) 



= 87rGc 2 /3o 



(11) 



where Rq denotes the scalar curvature corresponding to 
the unperturbed metric, /o = f(Ro), f R = df(Ro)/dRo 
and prime means derivative with respect to time rj. A 
very useful equation to use in the following calculations 
is the (fTTj) — (fTU)) combination 

,2\„2 



2(1 + f R ){-H' + U 2 ) + 2Hf' R - = 87rGp (l + c 2 s )a 



(12) 



Finally we have the conservation equation: 

P ' + i(l + c 2 )Hp = (13) 

Using the perturbed metric @ and the perturbed 
energy-momentum tensor the first order perturbed 
equations, assuming that the background equations hold, 
may be written as: 

(1 + f R )6GZ + (Rot + V^V, - 6£n)f RR 6R + 
[(5^ a )V,V a - S^5g^)W a W ]f R - 
Ig^iSTlu) -K9t \ST} a )]&r f R = -(toGST? (14) 

where f RR = d 2 /(i? )/di? 2 , □ = V Q V Q and V is 
the usual covariant derivative with respect to the un- 
perturbed FRW metric (see [l6| for perturbed metric, 
connection symbols and other useful perturbed quanti- 
ties). Notice that unlike the ordinary Einstein- Hilbcrt 



The Riemann tensor definition is ' va ^ = dgT{i a — daT^ 



in ie 
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case, with second order equations, this is a set of fourth- 
order differential equations. By computing the covariant 
derivative with respect to the perturbed metric V of the 
perturbed energy-momentum tensor Tjf, we find the con- 
servation equations: 



Let us consider equations (|16)) and (jl8|) for a matter 
dominated universe i.e. c 2 = 0, and combine them to 
express the potentials $ and \& in terms of {$', 8, 8'} 
by means of algebraic manipulations. The resulting ex- 
pressions are the following 



= 



(15) 



which do not depend on f(R). 

For the linearized Einstein's equations, the compo- 
nents (00), (ii), (Oi) = (iO) and (ij), where i,j — 1,2,3, 
i =/= j , in Fourier space, read respectively: 

+ [ -A: 2 ($ + *)-3H($' + *') + (3H / -6W 2 )$- 

m'^ ] +f R (-m<f> + m^> - W) = 2ps (16) 

(1 + /h) [ + + 3W($' + *') + + (W + 
2Ti 2 )^ ] +/^(3W$ - + 3$') + /r(3$ - *) = 
2c 2 JS (17) 

(l + /i ? )[$ , + *' + H($ + *)] + /^(2$-*) = 

-2p(l + (18) 



- p(7b){[ 3 ( 1+ ^w +$ ') 

+ + 2/551(1 + - + [(1 + + *') 

^(<5' - 3*')][(1 + /«)(-^ 2 - 3W0 + 3/^]} (24) 



and 
# = 



$ _ m = — -l^-SR 



l + /fi 



(19) 



where i5i? is given by 
2 



<5fl = - [3*" + 6(W + H 2 )$ + 3«($' + 3*') ~ 



fc 2 ($-2*)] 



(20) 



Finally, from the energy-momentum tensor conserva- 
tion (fl"5]) , we get to first order: 



3*'(1 + c 2 ) -5' + k 2 {l + c 2 s )v = 



(21) 



and 



" s <5W + ftu(l - 3c z s ) = (22) 



1 + c\ 



for the temporal and spatial components respectively. 

In a dust matter dominated universe, i.e. c 2 = 0, ^1 
and (|22p can be combined to give 



5" + ft<S' + k 2 <f> - 3*" - m^' = (23) 
which will be very useful in future calculations. 

EVOLUTION OF DENSITY PERTURBATIONS 

Our pourpose is to derive a fourth order differential 
equation for matter density perturbation 8 alone. This 
can be performed by means of the following process: 



P^-^y {[-3(1 + + *') - 3/^*' 

- 2p8][{l + f R )H + 2f R ] - [(1 + /«)($' + vl/') 

(25) 

where 

V(H,k) = -6(l + f R ) 2 H 3 + m[f'£ + 2{l + f R ) 2 H'] + 
3(l + fR)f' R (-2H 2 + k 2 +H') (26) 



The second step will be to derive equations (|24|) and 
with respect to ij and obtain and ^' algebraically in 
terms of {$", <5, 5', <5"}. These last results can be sub- 
stituted in equations p^|) and (fT5|) to obtain potentials $ 
and * just in terms of {$", <f", <5, <5', 5"}. So at this stage 
we are able to express, but we do not do here explicitly, 
the following 

$ = 8, 8', 8") 

* = 9 8, 8', 8") 

*' = *'($", <5', <5") (27) 

where we mean that the functions on the l.h.s. are al- 
gebraically dependent on the functions inside the paren- 
thesis on the r.h.s. 

The natural reasoning at this point would be to try 
to obtain the potentials second derivatives in 
terms of {6,6', 6"} by an algebraic process. The chosen 
equations to do so will be (|23[) and (|19p first derivative 
with respect to rj. In (f2"3"|) it is necessary to substitute $ 
and by the expressions obtained in (|2~7| whereas (fTS"]) 
first derivative may be sketched as follows 



fllR 

1 + fn 



6R' + 



[rrJr - [rrQ- + fn) 



(1 



6R 



Jr) 2 

(28) 

Before deriving, we are going to substitute \P" 
that appears on (fi"9|) by lower derivatives potentials 
{$, vp, <]/'}, 8 and its derivatives. To do so we consider 
(fTS)) and (fT5|) first derivatives with respect to rj where the 
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quantity v has been previously substituted by its expres- 
sion in (j2Tj) . Following this process we may express 
as follows 



if" = 6, 5', 5" 



(29) 



and now substituting in (|19p we can derive that equation 
with respect to rj. Solving a two algebraic equations sys- 
tem with equations (f23f and (|28]) and introducing (|27|) 
we are able to express {$", in terms of {5, <5', 6", 6"'}. 

$" = ; = *"(5,6',6",6'") (30) 



We substitute the results obtained in (|30|) straight- 
forwardly in (|2"T)) in order to express {$, <I>',\E''} in 
terms of {8, 8' , 5", 5'"}. With the two potentials and its 
first derivatives as algebraic functions of {6,S',6",S"'}, 
we performe the last step: We consider $(5, 8, 8" , 8"') and 
derive it with respect to rj. The result should be equal to 
$'(5, 8, 8", 8'") so we only need to express together these 
two results obtaining a fourth order differential equation 
for 8. Note that this procedure is completely general to 
first order for scalar perturbations in the metric formal- 
ism for f(R) gravities. 

Once this fourth order differential equation has been 
solved we may go backwards and by using the results for 8 
we obtain {<£>", from (|30j) as functions of time. Anal- 
ogously from (|27|) the behavior of the potentials {$, 
and their first derivatives could be determined. 

The resulting equation for 8 can be written as follows: 

faj8 m + fcjS'" + (a 2 , EH + 02,f)6" + Keh + Pij)6' + 
(«o,eh + Poj)8 = (31) 

where the coefficients fiij (i = 1, ...,4) involve terms 
with f' R and i.e. terms disappearing if we take fn 
constant. Equivalently, oj^eh (i — 0,1,2) contain terms 
coming from the linear part of /o in Rq . 

It is very useful to define the parameter e = TL/k since 
it will allow us to perform a perturbative expansion of the 
previous coefficients a's and /3's in the sub-Hubble limit. 
Other dimensionless parameters which will be used are 
the following: n t = /U l+1 (i = 1,2,3) and ft = 

fT/WfR) (7 = 1.2). 

Expressing the a's and /3's coefficients with those di- 
mensionless quantities we may write 



where two consecutives terms in each serie differ in e 2 
factor. The expressions for the coefficients are too long 
to be written explicitly. Instead, in the following sections 
we will show different approximated formulae useful in 
certain limits. 



VIABLE f(R) THEORIES 

Results obtained so far are valid for any f(R) theory. 
However, as mentioned in the introduction, this kind of 
models are severely constrained in order to provide con- 
sistent theories of gravity. In this section we review the 
main conditions Q: 



1- Irr > for high curvatures [17J. This is the re- 
quirement for a classically stable high-curvature regime 
and the existence of a matter dominated phase in the 
cosmological evolution. 

2. 1 + fn > for all i?o- This condition ensures the 
effective Newton's constant to be positive at all times 
and the graviton energy to be positive. 

3- /r < ensures ordinary General Relativity be- 
haviour is recovered at early times. Together with the 
condition fjm > 0, it implies that fn should be negative 
and monotonically growing function of i?o in the range 
-K f R < 0. 

4- |/r| <C 1 at recent epochs. This is imposed by 
local gravity tests [I3| , although it is still not clear what 
is the actual limit on this parameter. This condition 
also implies that the cosmological evolution at late times 
resembles that of ACDM. In any case, this constraint is 
not required if we are only interested in building models 
for cosmic acceleration. 



EVOLUTION OF SUB-HUBBLE MODES AND 
THE QUASI-STATIC APPROXIMATION 



Q?i, EH 

0u 



E a ?EH i = 0,l,2 
i=i 

3=1 
8 

E^3 ^0,1,2 

3=1 



(32) 



We are interested in the possible effects on the growth 
of density perturbations once they enter the Hubble ra- 
dius in the matter dominated era. In the sub-Hubble 
limit e <^ 1, it can be seen that the P^j and co- 
efficients are supressed by e 2 with respect to fh,fi Pij 
and /3o,/> i-e-, m this limit the equation for perturbations 
reduces to the following second order expression: 



5 



s" + m 1 



(1 + fiifn 2 {-\ + ki)(2ki - k 2 ) - |&/^(k 2 - 2)fc 8 87rGp a 



(l + /i0 5 (-l + Ki) + H/ 4 



(l + /ii)(K2-2)fe 8 



(33) 



r 



where we have taken only the leading terms in the e ex- 
pansion for the a and p coefficients. 

This expression can be compared with that usually 
considered in literature, obtained after performing strong 
simplifications in the perturbed equations - (|16p . (|17|) . 
(HHJ), ©, (EH) and d22J) - by neglecting time derivatives 
of $ and * potentials, (see [H|), Thus in [nj and [H 
they obtain: 



(5 



1 + 4 



Irr 



a 1 1+/h 



P<5 



1 , o/c 2 /flfi 1 I f 



= 



(34) 



This approximation has been considered as too aggres- 
sive in [111 ] since neglecting time derivatives can remove 
important information about the evolution. 

Note also that there exists a difference in a power k 8 
between those terms coming from the / part and those 
coming from the EH part in (|33|) . This result differs from 
that in the quasi-static approximation where difference is 
in a power k 2 according to (|3"4"1) . 
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Figure 1: 8 k with k = 0.2hMpc _1 for ftest (R) model and ACDM. 
Both, standard quasi-static evolution and equation J33D have been 
plotted in the redshift range from 100 to 0. 



In order to compare the evolution for both equa- 
tions, we have considered a specific function ftest{R) — 
— 4i? ' 63 , where H$ units have been used, which gives 
rise to a matter era followed by a late time accelerated 
phase with the correct deceleration parameter today. Ini- 
tial conditions in the matter era were given at redshift 
z — 485 where the EH part was dominant. Results, for 



k = 600H are presented in figure (1). We see that, as ex- 
pected, both expressions give rise to the same evolutions 
at early times (large redshifts) where they also agree with 
the standard ACDM evolution. However, at late times 
the quasi-static approximation fails to correctly describe 
the evolution of perturbations. 

Notice that the model example satisfies all the viability 
conditions described in the previous section except for the 
local gravity tests. As we will show in the next section, it 
is precisely this last condition \fn\ <C 1 what will ensure 
the validity of the quasi-static approximation. 



Recovering the quasi-static limit 

We will now restrict ourselves to models satisfying all 
the viability conditions, including \fn\ <C 1. 

In Appendix I we have reproduced all the a's and the 
first four /3's coefficients for each 5 term in (j3Tj) . These 
are the dominant ones for sub-Hubble modes (i.e. e<l) 
once the condition |/^| <C 1 has been imposed. Thus, 



keeping only Ej=i ^=o,...,4,/ and 01 
evant contributions for the general coefficients, the full 
differential equation (|31| can be simplified as 



.X 1 ) 



as the rel- 



(35) 



c 4 5' lv + c 3 S + c 2 5 + aS + cqS = 



where coefficients c's are written in Appendix II. 

We see that indeed in the sub-Hubble limit the C4 and 
C3 coeficients are negligible and the equation can be re- 
duced to a second order expression. 

As a consistency check, we find that both, in a matter 
dominated universe and in ACDM all f3 coefficients van- 
ish identically since fi, fa = 0. For these cases, equation 
(f3Tj) becomes equation (j6]) as expected. For instance, in 
the pure matter dominated case, coefficients k's are con- 
stant and they take the following values Ki = —1/2, 
K2 = 1/2, k 3 = —3/4 and K4 = 3/2. 

Another important feature from our results is that, 
in general, without imposing |/#| <C 1, the quotient 
(a^EH + Pij)/ ( a 2,EH + @2,f) is not always equal to H. 

In fact only the quotients Q^eh / ^2^eh and Pi^flP^f are 
identically equal to TL what is in agreement with the 8' 
coefficient in (|6]). However for our approximated expres- 
sions it is true that c\/c% = H. 

From expressions in Appendix II, the second order 
equation for 5 becomes 



6 



d" + H6'-- v , V ^l-/viH 2 J = 36 

' Gfnnk 2 , 5/1 /-i 24 M f S/aRfc 2 , 5 fi , /i 24 -l+«i \1 
~~ ^ 1" Si 1 ~ V 25 =2+7^ J L ^ h al 1 + Y 1 ~ 25 -2+k 2 )J 



3 

which can also be written as 



. ( 6f RR k 2 , 9\2 _ 81 i 9 2ki-k 2 

^ + ^-- V a2 - 4j 2 16 2 - 2+K2 (l- Kl )H 2 5 = (37) 

I a 3 V 2> _ X + b ^2+7^ 



Note that the quasi-static expression (|3"4")l is only recov- 
ered in the matter era (i.e. for H = 2/rf) or for a pure 
ACDM evolution for the background dynamics. Never- 
theless in the considered limit | fn |<C 1 it can be proven 
using the background equations of motion that 

1 + Ki - k 2 « (38) 

and therefore 2k\ — K2 ~ — 2+K2 ~ — what allows to 
simplify expression (|37|) to approximately become (f34|) . 
This is nothing but the fact that for viable models the 
background evolution resembles that of ACDM [j| . 

In other words, although for general f(R) functions 
the quasi-static approximation is not justified, for those 
viable functions describing the present phase of acceler- 
ated expansion and satisfying local gravity tests, it gives 
a correct description for the evolution of perturbations. 



SOME PROPOSED MODELS 

In order to check the results obtained in the previous 
section, we propose two particular f(R) theories which 
allow us to determine - at least numerically - all the quan- 
tities involved in the calculations and therefore to obtain 
solutions for (|3ip. As commented before, for viable mod- 
els the background evolution resemble that of ACDM at 
low redshifts and that of a matter dominated universe at 
high redshifts, i.e. the quantity (R + f(R))/R tends to 
one in the high curvature regime. Nevertheless the f(R) 
contribution gives the dominant contribution to the grav- 
itational action for small curvatures and therefore it may 
explain the cosmological acceleration. For the sake of 
concreteness we will fix the model parameters imposing 
a deceleration parameter today qo ~ —0.6. 

Thus, our first model (A) will be: f(R) = c x RP . Ac- 
cording to the results presented in [l2| and viable 
models of this type include both matter dominated and 
late-time accelerated universe provided the parameters 
satisfy ci < and < p < 1. We have chosen c\ = —4.3 
and p = 0.01 in Hq units. This choice does verify all 
the viability conditions, including \fji\ <C 1 today. For 
the second model (B): f(R) = e ; we have chosen 



ci = 2.5 -10 4 , e\ = 0.3 and c-2 = —0.22 also in the 
same units. 

For each model, we compare our result (f3"6"]) with the 
standard ACDM and the quasi-static approximation (f34|) 
(see Figs. 2 and 3). In both cases, the initial conditions 
are given at redshift z = 1000 where 6 is assumed to 
behave as in a matter dominated universe, i.e. 5k(ri) cx 
a(r)) with no k-dependence. We see that for both models, 
the quasi-static approximation gives a correct description 
for the evolution which clearly deviates from the ACDM 
case. 

In figure (H]) the density contrast evaluated today was 
plotted as a function of k for both models. The growing 
dependence of 8 with respect to k is verified. This mod- 
ified fc-dependence with respect to the standard matter 
dominated universe could give rise to observable conse- 
quences in the matter power spectrum, as shown in 
and could be used to constrain or even discard f(R) the- 
ories for cosmic acceleration. 




Figure 2: 8^ with k = 1.67hMpc _1 for f(R) model A evolving 
according to l|36|l . ACDM and quasi-static approximation given by 
equation (134 6 in the redshift range from 1000 to 0. The quasi- 
static evolution is indistinguishable from that coming from ((36} , 
but diverges from ACDM behaviour as z decreases. 
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Figure 3: 8 k with k = l.GThMpc" 1 for f(R) model B evolv- 
ing according to H36I I. ACDM and quasi-static evolution given by 
equation l|34| l in the redshift range from 1000 to 0. The quasi- 
static evolution is indistinguishable from that coming from (1361 . 
but diverges from ACDM behaviour as z decreases. 



\fa\ <C 1 today, then the perturbative calculation for 
sub-Hubble modes requires to take into account, not only 
the first terms, but also higher-order terms in e = TL/k. 
In that case, the resummation of such terms modifies the 
equation which can be seen to be equivalent to the quasi- 
static case but only if the universe expands as in a matter 
dominated phase or in a ACDM model. Finally, the fact 
that for models with |/#| <C 1 the background behaves 
today precisely as that of ACDM makes the quasi-static 
approximation correct in those cases. 
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Figure 4: Scale dependence of 5^ evaluated today (z = 0) for k/ Hq 
in the range from 1000 to 40000. 



CONCLUSIONS 

In this work we have studied the evolution of matter 
density perturbations in f{R) theories of gravity. We 
have presented a completely general procedure to obtain 
the exact fourth-order differential equation for the evo- 
lution of perturbations. We have shown that for sub- 
Hubble modes, the expression reduces to a second order 
equation. We have compared this result with that ob- 
tained within quasi-static approximation used in the lit- 
erature and found that for arbitrary f(R) functions, such 
an approximation is not justified. 

However, if we limit ourselves to theories for which 
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APPENDICES 



Appendix I: a's and 0's coefficients 



Coefficients for 5™ term: 



Coefficients for 5"' term: 



B {2) ~ 

fl(3) ~ 
P 3 ,/ - 

Pa,f — 

Coefficients for S" term: 



Coefficients for 6' term: 



0(1) 
8 {2) 



8/^(1 + fiiffU 2 
72/ 3 / 3 e 4 (-2 + k 2 ) 
216f 2 R f 2 e 6 (-2 + K 2 ) 2 
2f6/fl/ie 8 (-2 + K 2 ) 3 



(39) 



8/iHl + /fl) /fWe [3 + / fl (3 + fx)} 

6/lA 2 We 4 {8/ 2 (-2 + k 2 ) + 4/i[12ki + 9k 2 - 2(9 + k 3 )]} 

-72/|/!We 6 (-2 + «2)[-4/ 2 (-2 + k 2 ) + A(19 - 23 Kl - 10k 2 



4k 3 



-216/ fl We 8 (-2 + K2 ) 2 [-2/ 2 (-2 + k 2 ) + fi(7 - 11/si - 4k 2 + 2/s 3 )] 



(40) 



a (1) 
a (2) 

"2, EH 

(3) 
a 2,EH 

3 {1) 
"2,/ 

/3 (3) 
/3 (4) 

P 2,/ 



= 432(1 + / H ) 10 W 2 e 8 (-l + Kl )(-2 + K2 ) 3 

= 1296(1 + .fe) 10 W 2 e 10 (-l + ki) 2 (-2 + k 2 ) 3 

= 3888(1 + f R ) 10 H 2 e 12 (-l + Kl ) 2 (-2 + k 2 ) 3 

- 8/ 4 (1 + fnfffH 2 

~ 88f 3 R f 3 H 2 e 2 (-2 + K 2 ) 

~ 24/i/ 2 W 2 e 4 (-2 + K2 )(-28 + 2/si + 13/s 2 ) 

~ 72f R f 1 H 2 e 6 (-2 + k 2 ) 2 (-U + 4«i + 5k 2 ) 



(41) 



l.EH 



= 432(1 + .fe) 10 W 3 e 8 (-l + Kl )(-2 + « 2 ) 3 
= 2592(1 + / fl ) 10 W 3 e 10 (-l + «i) 2 (-2 + k 2 ) 3 
= -7776(1 + f R ) 10 H 3 e 12 (-l + mfi-2 + n 2 f 

~ sf R (i + .f R ) e f*n 3 

»;;. ~ 88/ 3 /jW(-2 + K2 ) 

24/£/ 2 H 3 e 4 (-2 + k 2 )(-28 + 2ki + 13k 2 ) 
72f R f 1 H 3 e 6 (-2 + k 2 ) 2 (-14 + 4«! + 5k 2 ) 



a (2) 

a l,EH 
a l,EH 

,(2) 



(42) 
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Coefficients for S term: 

"Seh = 432(1 + /fi) 10 H 4 e 8 (-l + Kl )(2 Kl -K 2 )(-2 + K2 ) 3 

ag H = 1296(1 + .f R ) w H i e 10 (-l + Kl ) 2 (-l + 4K 1 -K 2 )(-2 + K2 ) 3 

4 3 eh = 3888(1 + f R ) w H i e 12 (-l + K l ) 2 {2Ki-K 2 )(-2 + K 2 f 

~ -y/4(l + / fl ) 5 / 1 4 H 4 [2 + / ii (2 + 2/ 1 -/ 2 -2 Kl )-2« 1 ] 

/3g ~ 112/3/ 1 3 H 4 e 2 (-l + Kl )(-2 + « 2 ) 

/3<g ~ 48/^/ 1 2 W 4 e 4 (-l + K 1 )(-2 + K2 )(-16 + 2 Kl + 7K 2 ) 

^5 ~ 144/ fl / 1 H 4 e 6 (-l + K i)(-2 + K2 ) 2 (-6 + 4K 1 + K 2 ) (43) 
Appendix II: c's coefficients 



c 4 = -f R h[-fRhk 2 -'AH 2 (-2 + K2 )Y 

c 3 = - 3W 2 (-2 + « 2 )]{/l/ 1 3 fc 4 + 6/ 2 W 4 (-2 + « 2 ) 2 + AW 2 (-2 + n 2 )[2f R f 2 k 2 + 

m 2 (-7 + ll/si + 4k 2 - 2/s 3 )] + 2f R f 1 H 2 k 2 (-6 + 6«i + 3k 2 - k 3 )} 
c 2 = [-/fl/ifc 2 - 3W 2 (-2 + K 2 )] 2 [/iA 2 fc 4 + h] R j{H 2 k 2 {-2 + k 2 ) + 6W 4 (-1 + K!)(-2 + k 2 )] 
ci = Tic 2 

c = \h\-\ + K^l-fnhk 2 - ?M 2 (-2 + K 2 )] 2 [2f 2 f 2 k 4 + 9f R f 1 H 2 k 2 (-2 + k 2 ) + 9W 4 (2 Kl - « 2 )(-2 + k 2 )\ 

(44) 



